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Abstract 

We find a set of new exact solutions of a quantum harmonic oscillator, which describes 
some wave-packet trains with average energy being proportional to both the quantum level 
and classical energy of the oscillator. Center of the wave-packet trains may oscillate like a 
classical harmonic oscillator of frequency to. Width and highness of the trains may change 
simultaneously with frequency 2lo as an array of breathers. Under some perturbations the 
wave-packet trains could transit between the states of different quantum numbers. We 
demonstrate analytically and numerically that the wave-packet trains can be strictly fitted 
to the matter-wave soliton trains observed by Strecher et al. and reported in Nature 417, 
150(2002). When the wave-packets breathe with greater amplitudes, they show periodic 
collapse and revival of the matter- wave. 
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As an elementary equation of quantum mechanics the Schrodinger equation is a hnearly 
partial differential one of the second order with variable coefficients CP-jS|- We well known that 
this equation cannot be exactly solved yet for most physically interesting systems, except a 
few systems with separation of variables such as hydrogen atom, harmonic oscillator and rigid 
rotator The quantum states described by its solutions with inseparable space-time 

variables are very important but difficult to find. The coherent state of a harmonic oscillator is 
a nice example of such states [Hl-IHl- To seek new inseparable exact solutions of a Schrodinger 
equation and to physically realize them are our main motivations in this paper. 

The preparation and measurement of quantum states are very hard even impossible for 
some microscopic systems [5]- JI]. Compared to this the physical realization and detection of 
the macroscopic and mesoscopic quantum states may be easier sometimes J2]; The Bose- 
Einstein condensate (BEC) just supplies such a macroscopic quantum system which 
can be use to test the Schrodinger quantum mechanics. For example, a harmonically confined 
BEC can be identified as a perturbed quantum harmonic oscillator, when the interatomic 
interaction is weak enough ITTj, ^H]- The weak atom-atom interactions were due to the small 
atomic samples jTH], 120] and short s-wave scattering length \a\ [21], [22] that can be controlled 
by the Feshbach resonance |2S]; |21]- 

A quite interesting phenomenon was experimentally observed that the weakly interacting 
BEC appears solitonlike behavior |2I], [22] • Although this was approximately explained by 
using nonlinear interaction [2^], [2S]; Strecher and coworkers said that "non- interacting solitons 
■ ■ ■ would be expected to pass through one another" in their experiment. The solitonlike 
behavior is also explored for ideal BEC gas with a = in Khaykovich's experiment. Another 
important fact is the finding of collapse and revival of the macroscopic matter-wave packets 
|27j . |28j . which can occur for very weak interaction with atomic samples composed of a few 
thousand particles [TH] . In this paper we shall report a set of new exact solutions of a quantum 
harmonic oscillator and the corresponding macroscopic quantum level. By using them we 
demonstrate that, analytically and numerically, the wave-packet trains governed by these exact 
solutions can be strictly fitted to the matter- wave soliton trains found by Strecher et al.. This 
result means that Strecher's " non- interacting solitons" had been observed by themselves. On 
the other hand, under some particular initial and boundary conditions, these solutions exactly 
describe the well-known collapse and revival of a weakly interacting BEC. 
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We consider a BEC consisting of identical Bose atoms and being transferred into a cigar- 
shaped magneto-optical trap. Dynamics of the system is governed by the Gross-Pitaevskii 
equation (GPE) [TTI, [IHl 
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ih-^ = v' + [^m^'y + ^T^^liy' + z') + goW\iJ, (1) 

where ujr and uj^ are the transverse and axial frequencies respectively, the interaction intensity 
is related to the s-wave scattering length a, atomic mass m and number of atoms through qq = 
ATrNh^a/m for the normalized wave-function ^p. The norm is the probability density and 



N\^p\ the density of atomic number. Setting Ir = Jh/(muJr), h = y^/i'm^x) and writing E, 



-'kin 



and Eint as the kinetic energy and mean-field interaction energy of the BEC, the relationship 
Eint/ Ekin ~ N\a\/ {lllxY^'^ expresses the importance of the atom-atom interaction compared to 
the kinetic energy J7]. For low particle number jTHj, [20] or short s-wave scattering length \a\ 
|21j . we can treat the interaction term as a perturbation ^H] and obtain the leading order 
solution of Eq. (1) to obey the linear Schrddinger equation of a harmonic oscillator. Assuming 
the leading order wave-function is in the form of separation variable = il)n{x,t)il)y{y)'il)z{z) 
and its transverse factor is in the ground state of harmonic oscillator, we find a new exact 
solution of the linear Schrodinger equation (see Appendix A) 

= i'n{x,t)tpy{y)ij^{z) = Rn{x,y,z,t)exp[iQn{x,t)], 



1/2. 



Here bo is an arbitrary constant, -f^n(0 denotes the Hermitian polynomial of variable ^, Cq = 
9{t)p^{t) = ABujx sin(a — /?) is a conserved quantity of the classical harmonic oscillator equation 
ip = —uj^ip with the complex solution ip = p(t) exp[i^(t)], p{t) and 9(t) are the amplitude and 
phase of the complex oscillator 

p = \l cos^iujxt + a) + B"^ cos^iujxt + 6 = arctan — ( ^ ~^ f^) /gN 

^ Acos{LJxt + a) 

with A, B, a and /3 being arbitrary constants adjusted by the initial conditions of the classical 
harmonic oscillator. Obviously, the exact solution is not a energy eigenstate, but denotes a 
new kind of coherent states. In mathematical point of view, it is a complete solution with 
independent constants A, B, a, P and bo. By adjusting these constants, we can use the exact 
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solution to describe some different experimental results. It is easily to prove the exact solution 
(2) obeying the orthonormalization condition (see Appendix A). 

Amplitude Rn of the exact solution (2) describes the wave-packet trains consisting of n + 1 
packets. By using Eq. (3), from ^ = we have orbit of the center of wave-packet trains 

Xc — —p{t) cos9{t) = —Acos{u;xt + a), (4) 

Co Co 

which is proportional to real part of the complex solution if for a classical harmonic oscillator 
of unit mass with amplitude Abo/co, frequency cux and initial phase a. We call average energy 
of the system in state (2) as the macroscopic quantum level, which reads (see Appendix B) 

1 c 

En ^ + n)- + ujr + -C2, n = 0, 1, 2, (5) 
^2 ^ Co Co 

where constant Ci = {f? + c^/ + p^uoD /2 is another conserved quantity of the classical complex 
oscillator (3), 60C2/C0 — 6q74^ci/[cq(A^ -|- B'^)] is proportional to square of the amplitude of Eq. 
(4) or energy of the classical oscillator. Here and the following, we adopt the natural unit with 
m — h — ujx — ^ such that uj^ is normalized by uj^- For given constants, Eq. (5) exhibits that 
the energy only depends on the quantum number n. It is quite interesting that the average 
energy (5) is proportional to both the quantum level and the classical energies of the harmonic 
oscillators governed by Eqs. (3) and (4). 

The function p{t) included in ^ describes not only the total width of a wave-packet train 
but also the width of each packet. The average width of the packets is p{t) / i^co. Same function 
in radical of Eq. (2) governs highnesses of every packets. These and Eq. (3) infer that the 
widths and highnesses may change simultaneously with frequency 2uj. When the changes of 
the widths and highnesses are small, behavior of the wave-packet trains seemingly to be an 
array of solitons. And larger changes of the widths and highnesses show collapse and revival of 
the wave-packet trains, like multiple breathers. The normalization condition implies that the 
broader wave-packet train is associated with smaller mean highness and the narrower wave- 
packet train corresponds to larger mean highness. In the process of propagation and breathing, 
the wave-packet trains may spontaneously transit from state of higher average energy to that 
of lower one. Some perturbations also could cause transition between the states of different 
quantum numbers. All of these will be numerically illustrated as follows. 

Propagation of the wave-packet trains: pure bright soUton trains A pure soliton 
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always keeps its shape and energy in propagation, which can be described by the exact solution 
(2) through selecting p as a constant to fix the width and highness of each wave-packet. To 
do this we only require limiting the constants in Eq. (3) to A = B and f3 = a — tt/2 such 
that p = A. Adopting such selection to the second of Eq. (3) and to the formulas of constants 
Ci leads to the function 9{t) = uj^t + a and constants Cq = A^u^, Ci = oj^.A? , = oj^jl. 
In this case the ground state of Eq. (2) with = is just the common coherent state of 
a harmonic oscillator |E]-|H]- Inserting these constants into Eq. (5) yields the average energy 
En = {l/2+n)uJx+uJr+b'^/{2A'^). Final term of the energy is just the energy of classical harmonic 
oscillator (4), namely ^cuKboA/co)'^ = bl/{2A'^). Therefore, the average energy equates to sum 
of the quantum level and classical energy. We take the parameter set n = 10, a = 0, /? = 
—7r/2, A = B = p = co = l, &o = ~5, Ur = 4:0ux and adopt the units of time, space and 
probability density as u"^, Ix and from Eq. (2) to make the plot of on xoy plane 
for different times. In 1-3 lines of Fig. 1, we show the motions of eleven pure solitons from 
if: = 0, Xc(0) = — 5/a; to uJxt = 7r/2, Xc(vr/2) = 0, and to u^t = vr, Xc(vr) = bl^- In the 
propagation, the soliton train keeps its shape and distance between two solitons. 




Figure 1: The probability density of the wave-packet train as pure multi-solitons. Taking the 
parameter set n = 10, a = 0, [3 = — 7r/2, A = B = p = cq = 1, &o = ~5, uj^. = 40a;a; and 
normalizing the time, space and probability density in units uj~^, and from Eq. (2) we 
make the plot of B^ on xoy plane for different times. In the first line of Fig. 1 the initial profile 
of the soliton train is exhibited on left side of the trap. The soliton train propagates to center 
of the trap at uj^t = 7r/2 and to right side of the trap at uj: = vr, as in the second and third 
lines of Fig. 1. 
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Breathing of the wave-packet trains: collapse and revival For a very small constant 
Bq, Eq. (4) indicates that the amplitude of wave-packet oscillation may be very small. When 
the constant bo is taken as zero, center of the wave-packet trains is fixed to Xc = and their 
energy is reduce to En = (| + ^ + by Eq. (5). In order to show the collapse and revival, 
we let constant B be much greater than A, namely B = 1 and A = 0.01 such that the highness 
and width of the wave-packet trains oscillate with greater amplitudes. The other parameters 
are taken as n = 10, a = 0, P = — vr/2, cq = AB = 0.01, Ur/uJx = ll/^j = 40. Using the units 
of space-time coordinates and probability density of Fig. 1, we numerically draw the plots of 
vertical view for the wave-packet trains as Fig. 2. 




-50 -25 25 50 

X 



Figure 2: Collapse and revival of the wave-packet trains described by Eqs. (2) and (3) with 
parameters n = 10, bo = 0, A = co = 0.01, -6 = 1, a = 0, (3 = —n/2, Ur/uJx = ll/lf. = 40. 
The space-time coordinates and the density are normalized in same units with Fig. 1. Different 
colors represent different highnesses of the wave-packets as in right side of each line of Fig. 2. 
The initially high and narrow wave-packet train is shown in the first line of Fig. 2. As time 
increases to u^t = it/A, the wave-packets collapse to highness of 10"^ order, as in the second 
line of Fig. 2. To the half period Uxt = 7r/2, the highness is reduced to 10~^ order, as in final 
line of Fig. 2. 

In Fig. 2 we exhibit that center of the wave-packet train is fixed at Xc = 0, width and 
highness of each packet are greatly changed with period vr. The first line of Fig. 2 is corre- 
sponded to the initial wave-packet train, which is higher and narrower. The second line of Fig. 
2 shows the highnesses of packets have been reduced to 10^^ order at uj^t = 7r/4, and total 
width of the wave-packet train is simultaneously raised to 38{lx)- The highness and width of 
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the wave-packet train are collapsed to 10~^ order and 45 (/xO respectively, when time equates to 
the half period, as in third line of Fig. 2. In the next half period, revival of the wave-packet 
train will occur, through an inverse process of that described by Fig. 2. Fixing constant B, the 
forms of Eqs. (2) and (3) means that highnesses of the wave-packets are inversely proportional 
to constant A. When this constant is taken as very small, A ^ 0, highness of the wave-packet 
near a; = will tend to infinity, resulting in the intermittent implosions of the BEC j^l] , [^01 ■ 

Strecher's matter-wave soliton trains Generally, the wave-packet trains governed by 
Eq. (2) will propagate and breathe simultaneously, and spontaneously transit sometimes. We 
shall demonstrate that these behaviors can be strictly fit to Strecher's matter-wave soliton 
trains. In the experiment reported by Strecher and coworkers [2T|, the ^Li atomic BEC is 
employed to create the soliton trains, by using a Feshbach resonance to manipulate the sign 
and magnitude of the s-wave scattering length a. For a small value of \a\ and treating the 
atom-atom interaction as a perturbation proportional to a, the non-interacting soliton trains 
may be generated. 



3a 3b 3c 




Figure 3: The Strecher's matter-wave soliton trains on xoy plane from Eqs. (2) and (3) with 
the parameters n = 10, A = co = 0.4624, B = 1, a = 0, (3 = — vr/2, 6o = —17.437, ujr/ojx = 
ll/ll = 40, = 20Hz, /x = 21.22/im for (3a) t = 0, (3b) uj = n/2 and (3c) = tt. The 
space-time coordinates and the probability density are normalized in same units with Fig. 1. 
At t = 0, center of the soliton train is localized at Xc(0) = — 17.437(/x.) = — 370/im, average 
width of the solitons is p(0)/yco = 0.68(Zx), as in Fig. 3a. As time increasing to u^t = n/2, 
Fig. 3b shows that the center of the soliton train arrives at Xc(7r/2) = and the average width 
of the solitons is increased to p(7r/2)/^/co = 1.471(/3;). By Fig. 3c we display that the soliton 
train has moved to another end of the trap and its shape is changed back to initial case at 

U^t = TT. 
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It is specially interesting to investigate the formation of the soliton trains. Our theory shows 
that the non-interacting bright soliton trains could be generated for a ^ by the excitation 
from external fields. This means that the bright solitons described by Eq. (2) cannot be created 
in the case a > of Strecher's experiment. Only when the s-wave scattering length changes sign 
from positive to negative, the generating condition of the non-interacting solitons a ~ can 
be reached. Setting At as the interval between the time the end caps of Strecher's experiment 
are switched off to the time when a changes sign (a = 0), Strecher et al. found that number of 
the solitons increases linearly with At. At At = 0, namely the case that the end caps are 
switched off at a = 0, four non-interacting solitons were observed [21]. The larger At 
corresponds with stronger disturbatnce and the later can cause the system to higher excitation 
state with larger mean energy given in Eq. (5). The larger quantum number n is associated 
with more solitons. Because the known wave-packet of a quantum harmonic oscillator does not 
vary its width, Strecher et al. think of the soliton trains with variable width to be nonlinear 
multi-solitons. They also expected the non-interacting solitons being simultaneously released 
from different points in a harmonic potential. Given the exact solution (2) and macroscopic 
quantum level (5), the above analysis reveals that the soliton trains found by Strecher et al. 
just are the non-interacting solitons expected by themselves. 

Initial number of the solitons may be greater than ten in the experiment, oscillating fre- 
quency of the center of soliton trains is about 20Hz for the period T ^ 310ms and amplitude 
of that is 370/im. This frequency was identified as the axial one, Ux, in the previous ana- 
lytical work [25 . Their radial frequency is about 40 times the axial one. From Fig. 4 of 
Ref. [21] we estimate that the maximum and minimum widths of the soliton trains are about 
310/xm and 140;um. These experimental data give limitations to the parameters in Eqs. (2) 
and (3) as = cj^/40 = 20Hz, = ViOlr = 21.22/im, Abo/co = -17.437/^ = -370/im, 
p(0)/y^ = 6.8/^; = 144.30/im, p(7r/2)/yco = 14.71/^ = 312.15/im. Under these limitations, we 
choose the parameter set n = 10, A = cq = 0.4624, B = 1, a = 0, (3 = — vr/2, Bq = —17.437 
and use the units of space-time coordinates and probability density of Fig. 1 to make 3D plots 
of the soliton train for the time Uxt = 0, 11/2 and vr, respectively, as Fig. 3a, 3b and 3c. These 
plots display that the soliton train localized at ends of the trap possesses minimum width. 
When it move to center of the trap, its width becomes maximum. The change of the width was 
explained as repulsive interaction among the solitons in previous work In motion of 
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the soliton train, its highness has only small change such that effect of the collapse and revival 
cannot be observed. The solitons at ends of any soliton train are higher and thicker compared 
to the other ones. The thickness of each soliton is shown in Fig. 4, which is the vertical view 
of Fig. 3. This graph is very like the figure 4 of Strecher' article |2I], so the former could be 
a good fit to the latter. The differences of highness and thickness have not been accurately 
distinguished in the previous experiment. If this is done, we can fit it better by applying Eq. 
(2) and adjusting the constants A, B, a, [3 and ho. 




Figure 4: The vertical view of Fig. 3 with different line being corresponded to Fig. 3a, 3b and 
3c respectively. Comparison between this with the figure 4 of Strecher's article exhibits good 
agreement between them. 

Spontaneous transitions of the wave-packet trains In Strecher's experiment on the 
matter-wave solitons, the trains with missing solitons were frequently observed, and this is 
resided in loss of condensed atoms. However, "it is not clear whether this is because of a 
slow loss of atoms, or because of sudden loss of an individual soliton" [21]. According to 
transition theory of quantum-mechanical states, the non-interacting solitons could transit from 
higher-energy state to lower-energy state. The state function (2) and energy (5) imply that the 
lower-energy state describes less solitons. Therefore, even if the condensed atoms propagate 
without loss of number, some solitons may be suddenly lost, through spontaneous transitions 
of the macroscopic quantum states. Just after a transition, of course, the highness and width of 
each soliton should increase, since the atoms in lost solitons have entered the remainder solitons. 
The spontaneous transitions may be random and can be caused by some perturbations. We 
assume that under perturbation of the interatomic interaction the state of eleven solitons in 
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Figure 5: Transition of the matter- wave soliton trains. If the state of energy Eiq is disturbed 
by the interatomic interaction, the spontaneous transitions to the states of less average energy 
could occur. In Fig. 5 we show that the train consisting of eleven solitons is transformed into 
that of seven solitons for n = Q. The different lines correspond to the times Uxt = 2tt, 2.5tt 
and Sn respectively. 

Fig. 3 transits to the states with n = 6 at uj-J, = 2tt as in the first line of Fig. 5 and propagates 
to ujxt = 2.577 and u^t = Sir as in the second and third lines of Fig. 5. 

Further considering the first order correction to the wave-packet trains from weak interaction 
and investigating transformation of the interaction from weak to strong will be very interesting. 
Because of the existence of arbitrary constants A, B, a, f3, 6o cind periodic functions p(i(:), 0{t), 
by using the complete solution (2) we can control motions of the BEC wave-packets. The 
theoretical control could indicate the directions of experimental operations that is important 
for real application, say, making an atomic soliton laser based on the bright soliton trains. In 
addition, the exact solution (2) could play an important role in treating various harmonically 
confined systems. For example, a single Paul trapped ion interacting with a harmonic potential, 
the state \I'i of two wave-packets is similar to the Schrodinger's cat state jT2], [TH] . 

Acknowledgments This work was supported by the NNSF of China under Grant No. 
10275023 and the NLMRAMP of China under Grant No. T152103, and by the Hubei Provincial 
Key Laboratory of Gravitation and Quantum Physics of China. 
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Appendix A: Derivation of the Exact Solution 



Wc adopt the natural unit with m = h = uj^ = 1 and insert ^„ = ipnix,t)ipy{y)iljz{z) = 
{^/^:lr)~^ expl—^y"^ + z^) /{2l'^)]ipn{x,t), a = into Eq. (1), producing the one dimensional 
equation of a harmonic oscillator 



dt ~ 2 dx^ ^ 



V'n, (^1) 



where units of t and x are lo~^ and respectively, LOr is normalized by lo^ and appearance of 
the later is only formal, since its value has been fixed to 1. Let the solution of Eq. (Al) be in 
the form 

V'n = an{mn{i) exp[6(t)a; - c(t)a;2 - /=^(t)/2], i = e{t)x - fit) {A2) 

with a(t), c{t) being the complex functions of time and e(t), f{t) the real functions. 

Applying Eq. (A2) to Eq. (Al), we arrive at the equation 

0^ H OH 
e^-^^ + 2{be - if + iex - 2cex) 

+2 \i— - iff + — - c-ujr+ {ib - 2hc)x + {2c^ -id- -ujI)x'^] = 0. (^3) 
Noticing the Hermitian equation d'^Hn/d^'^ — 2^dHn/d^ + 2nHn = 0, Eq. (A3) implies 
id = 2c^ - ujI/2, ib = 26c, ie = 2ce - e^, 

i/ = be — e^f, idn/ttn = if f — b"^ /2 + c + cUr + ne^ . (^44) 

The first of Eq. (A4) is a complex Riccati equation, which can be changed into a complex 
equation of a classical harmonic oscillator 

'f = -^Iv, (^5) 

through the function transformation c = (p/{2i(p). The general solution of Eq. (A5) is well- 
known that 

if = A cos{u!xt + a) + iB cos{u!xt + /5) = p{x, t) exp[i6'(a;, t)] , (^46) 

where A, B, a and /3 are arbitrary constants, the real functions p{x,t) and 9{x,t) have been 
given by Eq. (3). Returning to the transformation between (p and c yields 

C = TTT- = -0 — I — . (A") 

2up 2 2p ^ ' 
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Substitution of Eq. (A6) into Eq. (A5) yields equations of the amplitude and phase as 



with the first integrations 



Co = p'Q = ABu^ sin(Q; - /?), Ci = {p' + c^/p^ + p'u;^)/2. 



Combining Eqs. (A7) and (A6) with Eq. (A4) and applying the relation (A9), we easily obtain 



cxi)(-/^?) 



cos^, 



A \ 
ttn — — exp ( — i ( — \-n)9 + LJrt — --^sin26' ). 



(AlO) 



^2 / 4co 

Inserting these into Eq. (A2) leads to the axial solution ipnix,t), and the normalization 



condition J {ipnl'^dx = Aly Cq^ J Hl{^) exp{—^'^)d^ = ttcq ^2"'n! = 1 gives the constant 
Aq = [^/(V7r2''n!)]V2 g^^h that Eq. (A2) becomes 

M^, t) = [ ^^ip^t^ ] "^H^^O exp [ - + ^Qn{x, t)] , {All) 

where the function Q{x, t) has been written in Eq. (2). Combining this with the wave-function 
of transverse dimension and letting the unit of spatial coordinate return to original one, x 
x/lxi we finally get the normalized 3D wave-function, as in Eq. (2). 



Appendix B: Proof of the Average Energy 



Employing the Dirac's symbols, ket and bra, from Eq. (A2) and the quantum-mechanical 
definition of average energy in state ■0„ we perform the calculation 



at a„ Hn ot, 

Noticing the orthonormalization condition {ipn\'4'n') = ^nn' and the formulas 

e = e{t)x fit), = y^V-n-i + 7(n + l)/2V'„+i, 

2eVn = ^n{n - l)iJn-2 + (2n + 1)V'„ + ^(n + l)(n + 2)V'n+2, 

1 dHfi _ Hn—l rz — 

Hn Hn 



{A12) 



12 



we continue to compute the average energy 

En = i{— - //) + i{ipn\V2ne^/e\ijjn-i) + i{il^n\hx - CX^IV'n) 



Applying Eqs. (A4) to the above equation results in 



En - (2n+l)c + ^,-^-(2c2-^)[/2+(i + n)]i + 26c^. 



(^13) 



Noticing Eqs. (A7), (A9) and (AlO) we have 



2 ^ Co 



C2 = - + f- - ^) cos^ e-^cose sin 



2 ^Co ^ P 

The third equation and Eqs. (A6) and (A9) imply 



C2 

e 



—p^ = - — cos^ 9 + ^(pcos^)^ — ^(psin^)pcos^ 

Co 2 Co Co 



Cl 



= cos -\ — 2 (P cos 6*) (psin6') — sin ^ — ^40;^ sin(a;j;t + o; 



Co 



Co 



Co . 



so that we get 



— [A^ cos^(a;^t + a) + B"^ cos^ {uj^t + /?)] 
Cq 

= — - + cos^(a;j;i + a) + —AB cos{uxt + (5) sin(a;-ri + a). 
2 Co Co 

Let the function of time be in the forms sin(2a;j;t) and cos{2ujxt), and identify the corresponding 
coefficients of both sides, producing 

A^Ci 

C2 



(A2 + 52)co' 

Combining these with Eq. (13) leads to 

K= (i + ^)- + ^r + -C2, n = 0, 1, 2, 



^2 ^ Co 

This is just the macroscopic quantum level (5). 
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